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The "ACFD-RPA" correlation energy functional has been widely applied to a variety of systems to success- 
fully predict energy differences, and less successfully predict absolute correlation energies. Here we present a 
parameter-free exchange-correlation kernel that systematically improves absolute correlation energies, while 
maintaining most of the good numerical properties that make the ACFD-RPA numerically tractable. The 
"RXH" kernel is constructed to approximate the true exchange kernel via a carefully weighted, easily com- 
putable radial averaging. Correlation energy errors of atoms with two to eighteen electrons show a thirteenfold 
improvement over the RPA and a threefold improvement over the related "PGG" kernel, for a mean absolute 
error of 13mHa or 5%. The average error is small compared to all but the most difficult to evaluate kernels, 
van der Waals Cq coefficients are less well predicted, but still show improvements on the RPA, especially for 
highly polarisable Li and Na. 
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Standard groundstate DFT methods 1 under local den- 
sity functional approaches such as the LDA, GGA or ex- 
act exchange (EXX) 2 approach have long been a fun- 
damental part of computational physics and chemistry. 
However, these functionals fail to reproduce important 
van der Waals (vdW) dispersion physics. 3 With increas- 
ing computer resources, the ACFD-RPA functional ap- 
proach, recently reviewed in Eshuis, Bates and Furche 4 
(EBF) and Ren et al, 5 has proved to be a versatile and 
fairly efficient 6-12 means of evaluating correlation ener- 
gies, albeit one that typically overpredicts absolute cor- 
relation energies due to its neglect of physics found in the 
exchange-correlation kernel / xc . 

The ACFD-RPA is typically employed as a "post 
groundstate" method, where a groundstate Kohn-Sham 
(KS) potential 1 is generated self-consistently in a stan- 
dard approach (eg. the LDA or EXX) and the orbitals 
are then used to refine the exchange-correlation energy 
without further self-consistency. Here the bare response 
function Xofa, x '] iw) is evaluated e.g. from the occupied 
and unoccupied orbitals of a groundstate KS system and 
used to evaluate the interacting response via 



Xa =Xo + Xo * w\ * xx 



(1) 



where * represents an integral/sum over x = ra. The 
effective interaction term w\ takes the form 

W X {x,x;iu>) =v^(\r ~r'\) + fl c (x,x'-iuj) (2) 

where vf (R) = X/R is a Coulomb potential and /^ c is 
the exchange-correlation kernel. Often f^ c (x , x' ; iu) is 
replaced by its zero frequency expression 



fl\x,x'-Q) = 



5V£ c (x) _ 5 2 El c 



5n(x') Sn(x)Sn(x) 



(3) 



or set to zero under the random-phase approximation 
(RPA). Finally the correlation energy can be evaluated 



through the so-called "ACFD" functional (see Section 2.1 
of EBF 4 for details) 
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Tr[Ca-Xo)*<f]. (4) 



This is then combined with the KS kinetic energy T KS , 
external energy E Ext and Hartree and exchange energy 
_E Hx to produce the electronic groundstate energy 



E n =T KS + E Ext + E Hx + E 



(5) 



By neglecting the exchange-correlation kernel / xc , the 
RPA has some considerable numerical advantages that 
are often exploited [see eg. Section 3 of EBF, 4 Ref. 13 and 
Ref. 5] to make calculations numerically efficient. How- 
ever, while the ACFD-RPA often predicts good energy 
differences, it is quite poor at predicting absolute corre- 
lation energies, especially in small atoms and molecules 14 
where it can overpredict correlation energies by as much 
as 150%. These inaccuracies come from failures in the 
short-range physics of interactions, arising from neglect 
of the / xc kernel. For example the opposite-spin pair- 
density n2CT^cr'(r, r') arising from the RPA can be neg- 
ative for r ~ r', a physically implausbile result that in 
turn contributes to overestimation of the absolute cor- 
relation energy. When taking energy differences these 
short-range errors tend to cancel each other out, which 
explains the RPA's success for these calculations. 

Various approximations to / xc , and other means to 
correct RPA absolute energies have been proposed and 
can be broadly divided into additive methods 14,15 that 
do not much affect dispersion interactions 3 and ker- 
nel methods 13 ' 16-24 that do. The most successful (and 
most numerically difficult) for atoms and molecules is 
the tdEXX 21,22 functional. Here the exact exchange ker- 
nel / x is used to produce very accurate correlation en- 
ergies, suggesting that the exchange kernel is sufficient 
for all but the most intractable problems. Indeed many 
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techniques 13,14 ' 19 22 base their corrections on exchange 
alone. 

Much of the numerical advantage of the RPA comes 
from the fact that w\ is: i) linear in A, ii) frequency 
independent, and iii) a function of \r — r'\ only. This 
allows one to perform the integration over A analytically, 
the integration over frequency efficiently, and, perhaps 
most importantly, reduce the complexity of convolutions 
(*) involving w\ by transforming to reciprocal space. For 
kernels with these three properties we can write the corre- 
lation energy as (in Hartree atomic units, used through- 
out this manuscript) 



g of the form 

g(x,x') = gaa'(R) = 



wfl 2 + (fc CTCT ,i?) 4 



(7) 



where R = \r — r'\ and the two parameters control the on- 
top curvature (c) and the range of the kernel (k). Clearly 
g(R — > 0) = c aa 'R 2 and g(R — > oo) = 1 as desired 
(except for H and He where c = k = correctly makes 
g oa = everywhere). The total interaction w\ = + 
f™ takes the desirable form 



w\(x, x 1 ) =A 



9*A\r-r'\) 



(8) 
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°° ^Tr 
o K 



{log(l + xo^i) - xowi } 

Wl 



(6) 



where xo(q,q'',iu) is a two-point function in q and 
q' while w x (<?,<?') = wi(q)S(q - q') and (<?,<?') = 

^jfq\ 5{q — q') are diagonal one-point functions. Clearly 
Wi(q) = vi(q) — |f in the RPA, but provided /f c is 
a function of \r — r'\ only, both iii and vi/ibi are di- 
agonal and (6) can be efficiently evaluated. Most ap- 
proximations to the true kernel cannot be written in this 
form, making routine application of kernel physics dif- 
ficult. The "Radial Exchange Hole" (RXH) kernel we 
shall develop in this manuscript aims to maintain these 
good numerical properties, while reproducing some good 
exchange kernel physics. 

An exchange kernel of intermediate accuracy and 
complexity is that proposed by Petersilka, Goss- 
mann and Gross (PGG). 20 They suggest approximat- 
ing the exchange kernel / x via the approach used 
by Krieger, Li and Iafrate 25 for the optimised ex- 
change potential. One can quickly 'derive' their ker- 
nel by first writing the Hartree and exchange en- 
ergy as E Hx = f ^0^g Hx (x,x')n(x)n(x r ), where g H x 
is the Hartree-exchange (Hx) pair-factor defined via 
the groundstate Hartree and exchange pair-density nux 
through gn x (x, x')n{x)n(x') = riYi x (x,x'). By as- 
suming 6gn x /5n m when evaluating ff c {x 1 x') = 

-\ one obtains the PGG kernel 



S 2 E Hx /[8n{x)6n(x')] - 
in the form f^ c (x,x') = A[gn x (x, x') — l]/|r — r'\ where 
.9x = [.9Hx — 1] < is the exchange hole pair-factor. This 
kernel has a number of good short-ranged properties com- 
ing from the use of the exchange hole g x . Notably it fully 
cancels the same spin Coulomb interaction for r = r' , 
and re- introduces it quadratically for small r — r'. For 
H and He the kernel reproduces the true f x = —8 aa /v c 
and spurious self-interaction is completely avoided. 

We now use an approximation similar to that of PGG 
to produce the more efficient RXH kernel. Here we aim 
to reproduce the good short-range physics of the PGG 
kernel, while replacing the two-point behaviour of (?Hx 
by a simpler dependence on the distance \r — r'\ only. 
We choose 26 a two parameter, "radial" model pair-factor 



We now introduce some "exchange hole" constraints to 
the model kernel. Rather than allowing the parameters 
c aa i and k aa i to be freely chosen, we instead ensure that 
the kernel has some integrated properties of the true pair- 
factor. Specifically we ensure that, for each spin pair 
ctct', the Hx pair number and groundstate Hx energy are 
reproduced via 

J drdr'{g(x, x')n(x)n{x') — nn x (x, x')] =0, (9) 



drdr' 



-\g(x,x')n(x)n(x') - n Hx (x,x')} =0. (10) 



The second constraint means that the 'derivation' of 
the PGG kernel given previously (where we assumed 
8g/Sn w in 5 2 E H */[8n8n]) can be made for the RXH, 
since £ Hx = / 



2i^n«'flHx = / ^% nn '9- The ker- 



nel might thus be considered a weighted radial averaging 
of the PGG kernel 

The integrals give the model pair-factor g some of the 
same physics as the true pair-factor g Hx - In particular 
they ensure that the dominant contribution from the ker- 
nel is short-ranged, with long-ranged effects only being 
felt in regions where the density (and thus response) is al- 
most zero. This in turn helps the RXH make its most im- 
portant changes to the correlation physics in regions with 
greatest contribution to the energy. Clearly the RXH ker- 
nel behaves quite differently to the the PGG kernel away 
from r = r'. For example, for electrons in different s 
shells (as in Li or Be) the PGG kernel gn x depends only 
on |r| and \r'\ and is exactly zero whenever |r| = \r'\, 
while the RXH g depends on \r — r'\ and is non-zero for 
electrons with |r| = |r'| but r ^ r' . However, if the effect 
of these differences on the correlation energy act in sim- 
ilar regions to those that dominate the exchange energy, 
the integral conditions (9) and (10) should ensure that 
the effect on the correlation energy is minimal. In fact, 
as results later show, the RXH kernel appears to better 
reproduce the true short-range physics of / xc than the 
PGG kernel, a somewhat surprising result. 

To test the RXH kernel, we apply it to atoms with two 
to eighteen electrons (He to Ar), using exact exchange 
groundstates calculations. Here we evaluate the ground- 
state properties under the LEXX 27 approximation to al- 
low for spherically symmeteric and spin homogeneous 
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densities. We use the same all-electron code previously 
employed to evaluate ISTLS correlation energies 23 with 
groundstate optimised effective potential approximated 
by the method of Kricgcr, Li and Iafrate 25 (KLI). The 
code is modified to reproduce the correct Hartree and 
exchange physics of ensembles via the method outlined 
in Ref. 27. We are thus able to reproduce groundstate 
properties, and in particular n 2 Hx, accurately for both 
closed and open shell systems. 

The linear-response calculation of xo involves utilising 
expansions on the spherical harmonics to produce cou- 
pled functions of radial coordinates r — \r\ and r' = \r'\ 
only. To generate xo we avoid sums over unoccupied 
orbitals and instead use Greens functions generated via 
a shooting method with errors coming only from the fi- 
nite radial grid with N r < 400 points. Energies were 
converged to under lmHa or within 0.5% (whichever is 
the larger). Given the simple radial geometry we do not 
transform integrals to reciprocal space, and the RPA, 
PGG and RXH methods all take the same O(N^) time. 

The parameters c aa i and k aa i of the RXH kernel (7) 
are evaluated by employing the secant method. Here, for 
a given k, we find c via (9), and k is then iterated to 
satisfy (10). Typically no more than ten iterations are 
required to accurately calculate k, although care must be 
taken to ensure convergence. Some example parameters 
(c aa /k aa ) are: Be) c = 0.127, k = 0.732, O) c = 0.667, 
k = 2.635, Na) c = 2.313, k = 3.475, Si) c = 4.583, 
k = 4.744, and Ar) c = 11.241, k = 5.692. 



TABLE I. Correlation energies in -mHa under the RPA, PGG 
and RXH approximations. Numbers in brackets are % errors. 
Groundstates are calculated in the LEXX approximation. 



System 


Exact 


RPA 


% 


PGG 


% 


RXH 


% 


He 


42 


84 


(100) 


45 


(7) 


45 


(7) 


Li 


45 


113 


(151) 


50 


(11) 


50 


(11) 


Be 


94 


181 


(93) 


104 


(11) 


106 


(12) 


B 


125 


224 


(79) 


112 


(-10) 


128 


(3) 


C 


156 


277 


(78) 


124 


(-21) 


155 


(-1) 


N 


188 


337 


(79) 


139 


(-26) 


175 


(-7) 


O 


258 


406 


(57) 


206 


(-20) 


243 


(-6) 


F 


324 


491 


(52) 


270 


(-17) 


304 


(-6) 


Ne 


390 


586 


(50) 


332 


(-15) 


369 


(-5) 


Na 


396 


611 


(54) 


332 


(-16) 


379 


(-4) 


Mg 


438 


674 


(54) 


375 


(-14) 


436 


(-0) 


Al 


470 


715 


(52) 


385 


(-18) 


467 


(-1) 


Si 


505 


772 


(53) 


401 


(-21) 


506 


(0) 


P 


540 


839 


(55) 


419 


(-22) 


552 


(2) 


s 


605 


903 


(49) 


474 


(-22) 


621 


(3) 


CI 


666 


982 


(47) 


527 


(-21) 


694 


(4) 


Ar 


722 


1072 


(48) 


579 


(-20) 


773 


(7) 


MAE(%) 




194 


(68) 


66 


(17) 


13 


(5) 



a From Ref. 28 



Absolute correlation energies — E c for the atoms are 



shown in Table I under the RPA, PGG and RXH ap- 
proximations. It is immediately clear that, with the ex- 
ception of Be, the RXH outperforms not only the RPA as 
expected, but also the PGG kernel. Here the mean ab- 
solute error of E c is reduced from 194mHa for the RPA 
and 66mHa for the PGG down to 13mHa for the RXH. 
Absolute percentage errors are reduced somewhat less, 
but even here RXH offers a threefold improvement from 
17% to 5% over the PGG kernel. For atoms with highest 
occupied orbitals in a p shell, the improvement is even 
greater still, with correlation energies approaching chem- 
ical accuracy. 

It is clear that the RXH must offer an unexpected im- 
provement over the PGG in the regions that dominate 
correlation energy calculations. It is likely that the spher- 
ical averaged form and integrated constraints correctly 
confine the kernel to be shorter-ranged than the PGG 
kernel. The dynamic physics will thus be reproduced 
better near the crucial "on-top" region. In the case of 
the transitions from p to s that dominate correlation en- 
ergies for outermost p shells the RXH kernel significantly 
outperforms the PGG kernel, adding further weight to 
the assumption that it is the short range physics that 
are most important. It would be very useful to com- 
pare the exact exchange kernel / x of Ne from a tdEXX 
calculation 21,22 with its RXH kernel. 



TABLE II. Same-species C& coefficients (Haajj) for spherical 
atoms in the RPA, PGG and RXH approximations. 





He 


Li 


Be 


Ne 


Na 


Mg 


Ar 


Exact" 


1.44 


1380 


219 


6.48 


1470 


630 


63.6 


RPA 


1.17 


500 


179 


4.98 


744 


481 


54.7 


PGG 


1.38 


1340 


277 


5.98 


1710 


743 


76.0 


RXH 


1.38 


1080 


231 


5.19 


1050 


487 


55.2 



a From lower bounds of Ref. 29 



We have so far demonstrated that the RXH kernel is 
good at reproducing the intra-species correlations that 
feed into the correlation energy. To test its ability to pre- 
dict inter-species correlations we also use the RXH kernel 
to evaluate van der Waals Cq coefficients for same-species 
pairs. These are notoriously difficult to reproduce accu- 
rately under standard DFT or ACFD type approaches. 
Here the correlation energy difference of two atoms at a 
great distance D obeys E vdw (D) = E C (D) - E c (oo) = 
— We calculate C% through the Casimir- Polder for- 
mula 

C 6 =^ J &loA(ilu) 2 , (11) 

where A(iu) = J dxdx' zz'xi(x, x') is the dipole polaris- 
ibility of a single atom. 

As shown in Table II the results are mixed, with the 
RXH kernel being only slightly more accurate than the 
RPA for the large Mg and Ar atoms, but otherwise out- 
performing it. Importantly it greatly improves on the 
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RPA for Li and Na, the systems with the greatest, and 
thus most dominant vdW C@ coefficients. We note that 
the evaluation of C§ coefficients is quite susceptible to 
even minor numerical issues and small variations in the 
groundstate orbitals, and these results may need to be 
better converged and tested for a comprehensive under- 
standing. 

In summary, we have produced a "Radial Exchange 
Hole" (RXH) exchange-correlation kernel that is both ef- 
ficient to calculate, and systematically improves on the 
RPA, at least for the atomic systems tested here. The 
kernel is derived from a similar approximation to that 
of PGG and is designed to model the exchange kernel 
/ x (r,r') via a simple dependence on \r — r'\ only. By 
writing the kernel in the form (7), constrained by (9) 
and (10) we are able to produce a kernel that is simple 
and parameter free. 

The RXH kernel offers improved correlation energies 
on all tested atoms (He to Ar) when compared with both 
the RPA and related PGG kernel (with the exception of 
Be). The mean absolute error of the atomic correlation 
energies is reduced to 13mHa (5%) from 66mHa (17%) 
and 194mHa (68%) respectively for the PGG and RPA 
kernels. Here the energy is dominated by the short-range 
physics correctly modelled by the RXH kernel. Results 
for vdW C% coefficients are more mixed, with predictions 
only slightly better than the RPA for Mg and Ar, but 
significantly better for Li, Be and Na. It is likely that 
longer-ranged kernel physics are involved in these calcu- 
lations, and these are more poorly predicted by the RXH 
approximation. 

While the RXH approximation derived and demon- 
strated here clearly performs well for correlation ener- 
gies of atoms, its immediate application to molecules 
and atoms may prove problematic. In particular the dis- 
sociation limit of different atoms will involve finding a 
pair factor g through (7) that is appropriate for both 
atoms. This problem could be remedied by treating the 
kernel as a sum over localised radial pair-factors [e.g. 

/ xc - E5©5(r)6 5 (r')^^pi with Q s (r i S) = 
for given regions S] , which we will explore in future work. 
Similarly in bulks, the form of the model (7) of the pair- 
factor is inappropriate, and an alternative approximation 
must be found. 

Overall, however, the RXH clearly represents an ap- 
proximation that is both numerically efficient, and more 
accurate than the RPA. Somewhat surprisingly it pro- 
duces more accurate correlation energies than the related 
PGG approximation, and many other approximations to 
/ xc . It is clear that constraining the kernel by its inte- 
grated properties allows it to accurately reproduce the 
important short-ranged physics that the RPA neglects, 
and this is likely to apply to other expressions for g. It 
should thus be possible to further refine the form of the 
kernel, perhaps via the addition of free parameters or 
a similar approach to include the correlation kernel, to 



improve both correlation energies and Cq coefficients. 
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